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Abstract 

We discuss the possibility to measure entropy of the system created 
in heavy ion colhsions using the Ma coincidence method. 



1 Introduction 

The assumption of thermodynamic equihbrium is one of the most commonly 
used when discussing the system created in central collisions of two relativis- 
tic nuclei. It is by no means obvious, however, that the equilibration actually 
can be achieved, since it is recognized as a process which may take longer 
time than the life time of the system in question. Be it or not, it is certainly 
important to verify if the created system is indeed in thermal equilibrium. 
To test this, one needs to check if the various measured quantities do satisfy 
the relations following from thermodynamics. In the present note we discuss 
the possibility of testing the relations pf] 
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and, with /i denoting the chemical potential, 
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which should be valid in any system at thermal equilibrium. 

Testing (|l|) requires measurement of the temperature T, the energy E, 
the number of particles n and the entropy S of the system in question. It is 
clear that by measuring the energies of the particles created in the collision 
we can measure the energy of the system. It is also generally accepted that 
by measuring the slope in the transverse momentum distribution we can 
measure the temperature^. The real difficulty is the measurement of entropy. 
In the present note we propose to adapt to this end the concidence method 
proposed some time ago by Ma . We also present Monte Carlo estimates of 
the feasibility of the method, based on a simple model. We conclude that the 
method has a large potential, as it requires much smaller number of events 
(~ y/no. of states), then the conventional approach. As a consequence its 
errors are significantly reduced compared to the simple-minded estimates. 
It is certainly worth to employ it in the present and future high energy 
experiments. 

2 Measurement of entropy by the coincidence 
method 

Below we give a summary of the idea presented in 0]. 

Ma proposes to count the pairs of configurations of the investigated sys- 
tem. Call Nc the number of pairs of "identical" configurations. Call Nt the 
total number of pairs of configurations. If all configurations considered are 
"equivalent" (i.e. if they correspond to the same conditions), then entropy 
is given by the formula 



The reason is that jf- is the volume in the phase-space occupied by the 
system. This can be seen as follows. 

^This requires correction to the effects of the hydrodynamic flow which seem to be 
under a reasonable control. 
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Suppose that the phase-space is divided in cells. Suppose furthermore 
that our system occupies F cells (with uniform probability). Each cell repre- 
sents a different state of the system (each cell has as many dimensions as is 
the number of variables describing the system). Our problem is to calculate 
T: S = log(r). Let us select randomly A/" configurations of the system (in 
general A/" ^ F, that's the main point). These configurations occupy some 
cells. The average occupation number of a cell is A/'/F ^ 1. Under this 
condition, the avergae number of pairs in the same cell is 



where Nf ^ N"^ is the total number of pairs selected. The total number of 
coincidences is the sum of (^ over all cells 

iv.^^(^)^ (5) 

hence 

and thus (D.| 

If the configurations are not equivalent, one has to divide them into 
classes: within each class they are now equivalent. If the probability dis- 
tribution of classes is -P(A), then 

The derivation is given in but it can be easily understood as a sum of 
the "average over classes" = J2xP{^) log (^^) ^^'^ "entropy of the 

distribution of classes" = — J2x ^og[P{X)]. 

The classes cannot be too small, so that number of configurations in each 
class is sufficient to make a reasonable statistics. 

This is what we retain from Ma. In the next section we present a sugges- 
tion how to apply this method to measure the entropy of a system of particles 
produced in high energy interactions. 

^The formula is only approximate. The exact formula is =^ ^^7^ which leads again 
to (1). 



3 



3 Application of coincidence method to mul- 
tiparticle production 

A natural possibility to apply the coincidence method to multiparticle pro- 
duction is to identify the configurations of the statistical system with the 
events observed in experiment. Once this is accepted, one can proceed as 
follows. 

(a) Select M events and split them into classes according to the total 
transverse energy E and multiplicity n recorded. The number of events in 
each class is denoted by M{E, n). 

(b) Define a "lattice" in momentum space (e.g. rapidity, azimuth and 
logpt, orlogEt) of individual particles Q 

Within each class: 

(i) Call the two configurations "identical" if they have the same occupa- 
tion numbers within the accuracy of the grid. The number of such pairs is 
denoted by Nc{E, n). 

(ii) Calculate the ratio (0), i.e.. 



where M{E, n) is the number of events in a given class. 

Actually, the condition that the events in one "equivalence class" must 
have strictly the same multiplicity (otherwise they could never be really iden- 
tical) could be relaxed, e.g. by accepting that in the definition of the "iden- 
tity" of the two configurations, the occupation numbers may differ by a small 
amount. 

If this procedure is going to have any sense, the results should depend on 
the lattice spacing a in a trivial way: when a — > a', S* ^ S" = 5* + log(a/a') 
(once the spacing is small enough). This must be checked, of course. 

Also question of size of the energy bins must be analysed. Again, the 
result should depend on the chosen size of the energy bins in a trivial way: 
if the size of the energy bin A changes into A' the entropy S changes into 
S' = S'-flog(A/A'). But this is delicate: the bins should be small enough so 
that within a bin the energy may be considered constant but large enough 

''Probably the better method is to transform the momenta into variables which give 
miiform distributions (see e.g. [^). 
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so that there is a reasonable statistics within each bin. Same apphes to the 
multiphcity. 

One sees from these arguments that in this way one can measure the 
entropy only up to an additive constant. Therefore the interesting thing is 
not to measure the absolute value of entropy but rather its dependence on 
energy or multiplicity. 

As we already noted in the Introduction, the measurement described by 
dl) allows one to perform a simple test of thermalization. When thermody- 
namics is valid, then the Eqs.(|l|,|) should be satisfied. Clearly the additive 
constant is irrelevant. One needs, however, a rather precise measurements be- 
cause otherwise the numerical estimates of the derivatives and ^^^J^'^^ 

oE on 

are not reliable. In the next section we show the results of a simple Monte 
Carlo estimate of the accuracies one can achieve. 



4 The classical gas of identical particles 

Of course the main problem is for how big systems the coincidence method 
works in practice. The number of states grows exponentially with the number 
of particles and the number of subdivisions. Therefore obviously there is a 
limit to what one can achieve with the finite computer. We will show here 
however that, the onset of the thermodynamic behaviour occurs still for the 
sizes where the coincidence method is feasible practically. 

Second, we will also show what must be the minimal size of the system 
for the continuum behaviour to set in. For very small number of subdivisions 
the coincidence method and our Monte Carlo are quickly convergent, however 
they converge to the correct description of the discrete problem, which is far 
from the continuum, hence not interesting. 

We consider the classical gas of noninteracing, nonrelativistic particles in 
d dimensions. Since, as mentioned earlier, the method may be applied to 
the transverse degrees of freedom only we prefer to retain the discussion in 
arbitrary dimensions. For the same reason we will use the number of degrees 
of freedom to characterise the size of a system. Of course N = d n in d 
dimensions 0. 

Since for noninteracting particles momentum and space degrees of free- 
dom factorize, we consider for simplicity only the momentum states. The 

^Provided only the momentum degrees of freedom are considered. 
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discretized expression for the number of states of degrees of freedom with 
the total energy E reads 



r{M,N)= Yl 1' (9) 

ni,...njv,n^+...n^=M 

where the momentum pi = arii with some discretization scale a. Accordingly 
2'mE = o?M , where the integer M labels the energy of the system and m 
denotes the mass of a particle. 
The generating function 

oo / oo \ ^ 

Zn{t) = E r''V{M,N) = E ^"M = exp(Anog c(r)), (10) 

Af=0 \n=-oo / 

factorizes and is expressed by a single sum c(r). The coefficients F can now be 
simply obtained by calculating recursively expansion of the /(r) = log c(r) 
from that of c(r), and subsequently expansion of Zj^{t) from that of /(t). 
This procedure provided us with the exact numbers for the density of states 
r(M, A^), which were used to benchmark the performance of our Monte Carlo. 
For large M and A^ the density of states reaches its continuum limit 

N/2AfiN/2-l) 

r{M,N)= , M,N large. (11) 

We will see later that this relation is rather well satisfied even for moderate 
values of M and A^. 

On the other hand, the thermodynamic limit, M, A^ oo, M/A^-fixed, is 
reached rather slowly. In this limit the entropy density Q scales depending 
only on the energy density e = M/N. 

1 1 

— logr(M, A^) ^ -(log (e) + log (27r) + 1), M, A^ ^ oo, e = M/N - const. 

^ (12) 
The purpose of this exercise is to see if the coincidence method can detect 
this behaviour. 



^Here and in the following, we will refer to the entropy per one degree of freedom as 
the entropy density. 
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4.1 Monte Carlo simulations 



In principle one should generate a sample of N configurations {ni, ?t,2, . . . , riArjfe, 
/c = 1, ...jA/" of, integer-valued, one-dimensional momenta ni, . . . ,nAr, which 
satisfy the energy conservation. For our purpose, however, the details of par- 
ticle kinematics, although practically cumbersome, are not relevant. In order 
to measure the coincidences, it is enough to label uniquely all multiparticle 
states and compare the labels. In this way the problem simplifies consider- 
ably, yet the essential question of the onset of the thermodynamic behaviour 
can be addressed. 

Consequently each Monte Carlo run consisted of a generation of a sample 
of N configurations, represented by integer indices, (Ji, J2, . . ., /a/"), 1 < /fc < 
r(M, A^), = 1, . . . , A/", uniformly distributed in the whole space of available 
states. Then we counted the number of coincidences Nc-, i.e., the number of 
pairs such that Ij = Ik- The estimate for the number of all states is 

then 

V = M{M -l)/Nc. (13) 

Moreover assuming the multinomial distribution of M integers among F bins 
we have calculated also the higher moments of the distribution of the number 
of coincidencies N^. In particular, the dispersion of reads Q 

a2[iV,] = 2<iV,>=^, (14) 
which gives for the relative error of the determination of F after M trials 



F M ' 

Therefore the estimate of the error, based on the MC data only, is 



^[F]/F = p/N,. (16) 

Eqs. (p!3|JT5|) show directly another advantage of the coincidence method. Namely, 
it works for much smaller number of trials (~ V^) than the standard ap- 
proach which measures average occupation of a single state. 

A sample of runs is summarized in Table 1. Exact results for F(M, A^) 
are also quoted. The last column gives the relative deviation of the current 

^After some approximations valid for 1 <C A/" <C F 
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estimate (col. 5) from the exact value. It should be compared with the 
estimate of the error based only on the Monte Carlo data, Eg. (|T6|) , given in 
column 6. The estimated error is steadily deacresing like l/M and actual 
deviation follows the suit albeit with some fluctuations. In all runs we have 
made (about 20 times more than shown in the Table) approximalety 30% of 
actual deviations were bigger that the MC estimate, as they should. Of course 
the formula (|15D is essential for planning future Monte Carlo simulations. It 
is interesting to note that the errors decrease as a number of trials and 
not as l/\/\Af. This is because the true random variable in this problem is 
the number of pairs, i.e. N"^. In particular the computing effort (counting 
pairs) grows like A/"^, and consequently the square root of the computational 
effort determines decrease of errors as it should. Altogether the Monte Carlo 
results are well under contol and show that the method is quite reliable. It 
is however practical only if the total number of states is less than several 
hundred milions. The last run shown in Table 1 lasted few hours on a 200 
MHz PC. This translates into A^, M <~ 25. We will discuss now if this is 
sufficient to see the onset of themodynamic properties. 



4.2 Results 

Figure 1 shows the entropy density as a function of a scaling variable e = 
M/N. Statistical errors of MC results (and the deviation from the exact 
discrete values given by r(M, A^)) are much smaller that the size of symbols. 
The data follow nicely the curves obtained from the classical formula in the 
continuum, Eq. (|Tl|) . Considered as a function of e and A^ they obviously 
show a substantial A^-dependence. The A^ varies from 8 (lowest curve) to 
24 in this plot. On the other hand, the deviation from the ultimate scal- 
ing limit, (Eq.([l^, the uppermost curve), is around 30% in the worst case 
(N=8,M=30). With A^ starting from 12, deviations from the infinite sys- 
tem are smaller than 20%. Note that A^ denotes the nubmer of degrees of 
freedom, which in d space dimensions corresponds to N/ d particles. 
As a second test we have checked a differential form of Eq. ( 



dlogV _ N 
dE 2E' ^ ' 

which, together with the equipartition of enegry, is the basis of the equilib- 
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rium thermodynamics. Changing the variable 2mE = a^M gives 
d\ogT _ dlogV dM _ dlogV M _ N 
dE ~ dM dE ~ dM 'E ~ 2E' 

or 



(18) 



dM 2M' ^ ' 

Finally after discretization of the derivative we obtain 

log (Ei^l±ll^) . (20) 

^ V r(M, N) J 2M+1 ^ ^ 

This equation is tested in Fig. 2, where a half of the inverse of the left hand 
side, as obtained from simulations, is ploted as a function of e. Solid line 
represents the right hand side [} Similarly to the previous case agreement is 
very good for > 12. It was neseccary to reduce MC errors to the level of 
l%-3% to achieve this agreement. Of couse this test is much more sensitive 
than the previous one since it requires precise measurement of the derivatives. 

To conclude, the coincidence method is satisfactory in practice for the 
number of degrees of freedom below ~ 25. This turns out to be sufficient 
to see the signatures of the thermal equilibrium. For more than 12 degrees 
of fredom the scaling of the entropy density is confirmed with the accuracy 
better than 20% . The saddle point relation dS/dE = 1/T is also very well 
reproduced. 

This work is supported in part by the Polish Committee for Scientific 
Research under the grants no. 2P03B 08614 and 2P03B 04412. 
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N 


M 


u 




r 


a/V 


5/V 






4 000 


218 


73 376. 


0.096 


0.140 






8 000 


1000 


63 299. 


0.045 


0.007 




6 


16 000 


3866 


66 214. 


0.023 


0.028 






32 000 


15884 


64 465. 


0.011 


0.001 




r 


64 416 






8 000 


30 


2 133 067. 


0.260 


0.100 






16 000 


124 


2 064 387. 


0.127 


0.065 


12 


12 


32 000 


516 


1 984 434. 


0.062 


0.024 






64 000 


2 110 


1 941 201. 


0.031 


0.001 






128 000 


8 358 


1 960 262. 


0.015 


0.011 




r 


1 938 336 






20 000 


4 


99 995 000. 


0.707 


0.615 






40 000 


16 


99 997 504. 


0.354 


0.615 




24 


80 000 


106 


60 376 604. 


0.137 


0.025 






160 000 


422 


60 663 128. 


0.069 


0.020 






320 000 


1596 


64 160 200. 


0.035 


0.036 




r 


61 903 776 






8 000 


8 


7 999 000. 


0.500 


0.077 






16 000 


30 


8 532 800. 


0.258 


0.015 




6 


32 000 


132 


7 757 334. 


0.123 


0.105 






64 000 


442 


9 266 824. 


0.067 


0.070 






128 000 


1842 


8 894 610. 


0.033 


0.027 


24 


r 


8 667 720 






40 000 


78 


99 997 504. 


0.354 


0.487 






80 000 


2100 


199 997 504. 


0.250 


0.025 




8 


160 000 


8334 


209 834 752. 


0.128 


0.076 






320 000 


33658 


200 783 680. 


0.063 


0.029 




r 


195 082 320 



Table 1: Monte Carlo results for r(M, A^) (col. 5) for different and M. The 
third and fourth column give the number of generated configurations A/", and 
the number of observed coincidences N^- In the last two columns we quote 
the Monte Carlo estimate of the relative error ,c.f. Eq.(|TB|), and the actual 
relative deviation 5/T = |r — r|/r from the exact value F also quoted in the 
Table. 
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Figure 1: Entropy density s = ^\ogW{M, N) vs. the energy density e = 
M/N . Black symbols represent our Monte Carlo results for N=8 (diamonds), 
12 (circles), 16 (boxes) and 24 (a triangle). Lower solid lines correspond to 
the continuum approximation, Eq.(|ll|), for each A^. The uppermost solid line 
represents the scaling, thermodynamical limit, Eq.(|T^. 



11 



-3L . . . . . . . . . ^ . . . . ^ . . . . J 
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Figure 2: Testing the relation (pOl). Half of the inverse of the finite difference 
(with respect to M) of the entropy logr(M, iV), as a function of e = M/N 
for N = 8 (diamonds), 12 (circles), 16(boxes) and 24 (a triangle). Solid line 
corresponds to the thermodynamical Imit. 
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